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ABSTRACT 

Simulations and analytic arguments suggest that the turbulence driven by magne- 
torotational instability (MRI) in accretion discs can amplify the toroidal (azimuthal) 
component of the magnetic field to a point at which magnetic pressure exceeds the 
combined gas + radiation pressure in the disc. Arguing from the recent analysis by 
Pessah and Psaltis, and other MRI results in the literature, we conjecture that the 
limiting field strength for a thin disc is such that the Alfven speed roughly equals 
the geometric mean of the Keplerian speed and the gas sound speed. We examine 
the properties of such magnetically-dominated discs, and show that they resolve a 
number of outstanding problems in accretion disc theory. The discs would be thicker 
than standard (Shakura-Sunyaev) discs at the same radius and accretion rate, and 
would tend to have higher colour temperatures. If they transport angular momentum 
according to an a prescription, they would be stable against the thermal and viscous 
instabilities that are found in standard disc models. In discs fuelling active galactic 
nuclei, magnetic pressure support could also alleviate the restriction on accretion rate 
imposed by disc self-gravity. 
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1 INTRODUCTION 

The standard theory of geometrically thin accretion discs 
(Shakura & Sunyaev 1973) has provided a successful frame- 
work for understanding the basic physical properties of these 
objects. And good progress toward understanding the main 
uncertainty in the theory, the viscosity, has been made since 
it was realized that the magneto-rotational instability (MRI) 
can drive dynamo action to maintain magnetic stresses (Bal- 
bus & Hawley 1998). However, several aspects of disc theory 
have remained problematic over the years. 

Four of these problems stand out. First, observed disc 
spectra in cataclysmic variables (CVs) deviate significantly 
from theoretical expectations, even when self-irradiation is 
taken into account. Moreover, different kinds of CV systems 
with similar accretion rates often show different spectra. Sec- 
ond, the discs in CVs and low-mass X-ray binaries (LMXBs) 
often appear to be geometrically thicker than expected, sug- 
gesting that some additional source of pressure supports the 
optically thick layers against gravity, over and above the 
usual mix of gas + radiation pressure. Third, disc models in 
which the viscous stress is roughly proportional to the total 
pressure ( "a- models" ) — an otherwise plausible assumption 
— are often predicted to be unstable to perturbations in 



temperature and surface density. Yet observations of many 
potentially unstable sources do not show clear evidence of 
such instabilities. Fourth, the accretion discs thought to fuel 
active galactic nuclei (AGN) are predicted to be gravitation- 
ally unstable on scales exceeding a fraction of a parsec, yet 
mass transport does not seem to be disrupted by fragmen- 
tation and star formation in many AGN. 

In this paper, we propose a mode of disc accretion in 
which the pressure of the organized toroidal magnetic field 
dominates over other forms of pressure. On the basis of pub- 
lished numerical simulations, we argue that the MRI-driven 
dynamo could amplify the toroidal field to a strength that 
greatly exceeds the gas pressure. Similar arguments were 
presented by Pariev, Blackman & Boldyrev (2003) and Blaes 
et al. (2006), but the former focused on the stochastic part 
of the field, and neither attempted to estimate the field 
strength from first principles. Arguing on the basis of anal- 
yses by Pessah & Psaltis (2005) and others, we suggest that 
MRI-driven turbulence could amplify the toroidal field to a 
maximum strength such that the associated Alfven speed is 
roughly the geometric mean between the Keplerian speed 
and the gas sound speed. Such a strong field can dominate 
the disc structure even in regions where radiation pressure 
dominates over gas pressure. We construct a model of a mag- 
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netically dominated a-disc according to this prescription, 
and show that the additional pressure support both thick- 
ens the disc and increases its colour temperature, stabilizes 
thermal and viscous instabilities, and ameliorates the effects 
of gravitational instability. 

The plan of the paper is as follows. In Section 2, we out- 
line the evidence for anomalous spectra and disc thicknesses 
in CVs and LMXBs. We defer brief discussions of the stabil- 
ity problems to Section 4. In Section 3, we discuss the numer- 
ical evidence supporting strong toroidal fields in discs, then 
focus on the Pessah & Psaltis (2005) analysis and results, as 
well as results in the literature for nonaxisymmetric MRI, in 
order to justify our estimate for the characteristic disc field. 
In Section 4 we adopt eq. ((4]) as the constitutive relation 
for the magnetic pressure, and construct a family of simple 
black hole disc models characterized by an a— viscosity un- 
der this assumption. We show that the resulting discs are 
hotter and thicker than standard discs. Most importantly, 
and in contrast to the inner regions of standard black hole 
disc models, they are not subject to thermal or viscous insta- 
bilities. In addition, in the outer regions of AGN discs, where 
self-gravity is important, they can carry much larger mass 
fluxes without fragmenting due to gravitational instability. 
Finally, we discuss our results and summarize our conclu- 
sions in Section 5. 



that these discs are optically thick. Thus, if the accretion en- 
ergy is deposited within the body of the disc (where most of 
the matter resides), computing the spectra should be closely 
akin to the well understood problem of computing the spec- 
tra of stellar atmospheres, especially given that in this case 
the temperatures and gravities (densities) involved are in 
the standard range (Wade & Hubeny 1998). 

However, it has long been known that this simple- 
minded approach does not give satisfactory results. Just 
from an empirical point of view it was already clear that 
there might be problem: in her review of IUE spectra, la 
Dous (1991) showed that in the range AA1200 - 3000 the 
continuum spectra of dwarf novae at the peak of an outburst 
are significantly "bluer" (i.e. have more flux at shorter wave- 
lengths) than the continuum spectra of non-magnetic nova- 
like stars. For the nova-like variables, which should most 
closely approximate steady-state discs, it turns out that it 
is difficult, if not impossible, to replicate the UV data us- 
ing steady-state disc models constructed from LTE stellar 
atmospheres (Wade 1988; Long et al. 1994; Orosz & Wade 
2003). In general, the models predict UV spectra that are too 
"blue," and in addition it is not possible to match both the 
observed flux and colour at the same time. Similar problems 
are found for dwarf novae at the peak of outburst (Knigge 
et al. 1997) 



2 PROBLEMS WITH CV AND LMXB DISC 
MODELS 

2.1 Spectral Properties 

In principle a steady-state, geometrically thin accretion disc 
should be relatively simple to understand. Since the accre- 
tion rate is independent of radius, R, the effective tempera- 
ture T c ff(R) is independent of the viscosity mechanism, and 
is given just in terms of system parameters (Pringle 1981) 
as 



r-m-T. (A)- [,.(-) 



-1/2' 



1 1/4 



(1) 



where R, is the radius of the inner disc edge, and the ref- 
erence temperature T„ is given in terms of the mass of the 
central object, M*, and the accretion rate, M, by 

where asB is the Stefan-Boltzmann constant. 

The astronomical objects to which this formula should 
be most straightforwardly applicable are the accretion discs 
in cataclysmic variables (CVs), and in particular the discs 
in nova-like variables, and dwarf novae in outburst. In the 
nova-like variables, and in the Z Cam variables in standstill, 
the discs are steady-state in that the timescale for varia- 
tion of the accretion rate (measured by system brightness) 
is much longer than the viscous timescale (as deduced em- 
pirically from the evolution timescale for dwarf nova out- 
bursts). Indeed, even in dwarf novae in outburst, and espe- 
cially in super-outburst, the discs are close enough to being 
steady state that their spectra should be readily calculable. 
Estimates of accretion rates and evolution timescales imply 



2.2 Disc thickness 

In the standard accretion disc picture, the thickness of the 
disc is governed by a balance between pressure and (a com- 
ponent of) gravity. To a first approximation, the scaleheight 
of the disc, H, is given by 



H ~ R 



(3) 



(Pringle 1981), where c s is the sound speed at the disc mid- 
plane and Vcf> is the azimuthal velocity. For a simple steady- 
state CV disc, the shape of the disc is slightly flared in that 
H/R oc R 1/8 (Shakura & Sunyaev 1973). 

In order to fit the eclipse profiles of discs in CVs, it is 
necessary to allow for the fact that the disc is not flat, and 
for the fact that the disc photosphere is usually at a height 
of a couple of disc scaleheights above the disc midplane. 
Eclipse fitting does allow some estimation of the disc thick- 
ness, usually in terms of the value of H/R at the outer disc 
edge. The observationally inferred values of disc thickness 
required to fit eclipse observations are typically a couple of 
times larger than the theoretical estimates (Robinson, Wood 
& Wade 1999; Shatter & Misselt 2006). 

There are other indications of material farther out of 
the plane than simple estimates would suggest, notably the 
"iron curtain" material seen both in quiescent discs (Home 
et al. 1994) and at high mass transfer rates (Baptista et 
al. 1998), and transient non-axisymmetric features, possi- 
bly spiral structure seen in tenuous layers of the atmosphere 
(Steeghs, Harlaftis & Home 1997), which may require an un- 
usually hot, or thick, component of the atmosphere and/or 
additional tidal heating at the disc edge (Ogilvie 2002). 

LMXBs (White & Mason 1985; Hakala, Muhli & Dubus 
1999) and the Algol binary W Cru (Pavlovski, Burki & Mim- 
ica 2006) provide more evidence that the outer disc edge 
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does not fit simply with standard disc theory, in that the 
outer disc edge seems to be thick and structured. In ad- 
dition, the apparent paucity of eclipses among LMXBs, or 
bulge sources, indicates that the disc manages to obscure 
the view of the central X-ray source from a larger fraction 
of the sky than might be expected (Milgrom 1978; Joss & 
Rappaport 1979). 

The disc rim may not be the ideal place to test standard 
disc theory, since the outer disc edge is subject to both tidal 
dissipation and the impact of the accretion stream (Lubow 
1989; Armitage & Livio 1998). Nevertheless, the soft X-ray 
transients do seem to indicate that in these objects the disc 
thickness is larger than expected at most radii, and not just 
at the outer disc edge. This is because, in order to explain the 
prolonged nature of the outbursts and the shape of the decay 
lightcurve, the simplest explanation appears to be that the 
disc is significantly irradiated from the centre (King 1998). 
According to standard disc models, such irradiation is not 
possible, since the inner parts of the disc screen the outer 
regions (Dubus et al. 1999). However, one can provide an 
adequate description of the outbursts of soft X-ray tran- 
sients by assuming that the disc is thicker than simple mod- 
els predict, allowing significant irradiation to occur (Dubus, 
Hameury & Lasota 2001). 



3 HOW STRONG ARE DISC MAGNETIC 
FIELDS? 

It seems likely that the problems discussed above (and sev- 
eral to be discussed later) could be resolved by considering 
the dynamical effects of magnetic fields in accretion discs. 
In standard disc theory, spectral predictions are based on 
two assumptions: first, that energy is deposited where most 
of the mass is, and second, that the density structure of 
the disc is smooth and homogeneous. However, it has long 
been realized that if the process that taps the shear energy, 
first converts it to magnetic energy, which is then dissipated, 
then much of the energy can be deposited in a small fraction 
of the mass (Lynden-Bell 1969; Meszaros, Meyer & Pringle 
1977). Moreover, it has also long been known that if a disc 
is strongly magnetic, it is likely to be highly inhomogeneous 
(Pustilnik & Shvartsman 1974). In recent models of spec- 
trum formation in the low/hard state of discs in LMXBs 
and for the hard X-ray flux in AGN, these ideas are taken 
to the limit in which all the accretion energy is assumed to be 
deposited in a low density, extended magnetic corona (see, 
for example, Zycki, Done & Smith, 2001; Done & Nayak- 
shin, 2001, Barrio, Done & Nayakshin, 2003; and references 
therein). 

Numerical simulations give some support to these ideas. 
Hirose, Krolik & Stone (2006) and Fromang & Nelson 
(2006), following on from earlier work by Miller & Stone 
(2000), present the results of shearing-box simulations of 
MRI-driven accretion discs, in which the vertical extent of 
the computational grid is large enough to encompass regions 
of low gas density. The typical structure that they find is a 
magnetic sandwich, with most of the mass forming a gas 
pressure-dominated disc in the plane, but with extensive, 
low density, magnetically supported layers on either side. In 
these simulations, most of the energy dissipation still occurs 
in the gas pressure-dominated central layer, and therefore 



at high optical depth. However, the optical depth of the ex- 
tended magnetic atmosphere is non-negligible, and this has 
two effects. First, the density at the effective photosphere 
is lowered, increasing the ratio of scattering to absorptive 
opacity, and so giving rise to a slightly harder spectrum 
(Blaes et al. 2006) . Second, the height of the apparent pho- 
tosphere is increased, making the disc appear thicker than 
a standard disc. 

In this paper we raise the possibility of there being a sec- 
ond equilibrium distribution for the magnetic flux in a disc, 
namely one in which the magnetic pressure (pb ~ B 2 /8ir) 
dominates at all heights in the disc (or at least through most 
of the pressure scale height) . The existence of such a distinct 
class of discs has been raised before by Shibata, Tajima & 
Matsumoto (1990). As we discuss below, such a possibil- 
ity would be difficult to simulate numerically because of the 
small lengthscales of the instabilities that maintain the mag- 
netic dynamo. It is also likely that the formation of such a 
strongly magnetized disc requires special conditions of some 
kind. For example, Pringle (1989) discusses the formation 
of a strongly magnetized region in an accretion disc bound- 
ary layer, where the strong shear leads to a much higher 
formation rate for toroidal flux. A disc might also become 
strongly magnetically dominated as a result of thermal in- 
stability, for example in the central regions of black hole 
accretion discs (Machida, Nakamura & Matsumoto 2006), 
or during the transit to quiescence in a cataclysmic variable 
disc (Tout & Pringle 1992). In these circumstances the disc 
is originally hot, with ps ~ ap g , where p g is the gas pressure 
and a is the Shakura-Sunyaev (1973) viscosity parameter. 
The disc then cools (reducing p but not pb) rapidly on a 
thermal timescale (~ 1/aQ, where Q is the angular speed). 
This would result in a disc in which /3 = p 9 /pb <S 1. 

Pariev et al. (2003) have considered the structure of a 
strongly magnetic disc in which /3<q~ 1 , but were unable 
to predict the strength of the field. In this paper, we attempt 
to quantify the characteristic strength of the toroidal mag- 
netic fields likely to dominate accretion discs, and to explore 
the consequences of these strong fields for vertical and ra- 
dial disc structure. We base our estimate of the limiting field 
strength on a recent analysis of the effects of strong toroidal 
fields on MRI by Pessah & Psaltis (2005), as well as other 
analyses in the literature. We assume that toroidal fields 
are amplified by dynamo action resulting from the turbu- 
lence driven by MRI, and are probably limited by buoyancy 
effects operating on some multiple of the dynamical (orbital) 
timescale. As long as MRI grows on a dynamical timescale, 
the azimuthal field, B$, continues to be amplified, but when 
the MRI growth rate is suppressed by magnetic tension the 
growth stops and saturates. 

MRI continues to operate on a dynamical timescale 
even when the toroidal magnetic pressure exceeds the gas 
pressure. Pessah & Psaltis (2005) show that the maximum 
growth rate for axisymmetric MRI modes begins to be 
severely suppressed only when the Alfven speed associated 
with the toroidal field, v\4> = (B^/Anp) 1 ^ 2 , exceeds the ge- 
ometric mean of the Keplerian speed vk and the gas sound 
speed c g . (Radiation pressure is not considered in their anal- 
ysis; we address this later.) When vai/> = (2c 3 vk) 1 ^ 2 , the 
growth rate of MRI vanishes. Although different types of 
axisymmetric and nonaxisymmetric modes may grow when 
the field strength is larger than this limit, we conjecture 
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the va</> cannot grow past this instability "bottleneck." We 
therefore adopt this limiting value for axisymmetric MRI as 
a measure of the characteristic magnetic pressure, at least 
when gas pressure dominates: 



Pb ~ pc 3 VK- 



(4) 



In this limit the magnetic pressure exceeds the gas pressure 
by a factor 



/r 1 



Pb 
Pa 



> l; 



(5) 



thus, magnetic pressure would dominate the structure of the 
disc. 



3.1 MRI and the Limiting Toroidal Field Strength 

In this section, we analyze the physics behind the limiting 
field strength, and show its relationship to other stability 
results in the literature. Pessah & Psaltis (2005) consider a 
differentially rotating, cylindrical equilibrium flow with no 
z-dependence. In their analysis, they make approximations 
that correspond to assuming that the flow is of uniform den- 
sity and contains a uniform poloidal field B z 7^ and a uni- 
form toroidal field B^, (their results can be generalized to in- 
clude a radial gradient of B$). The angular velocity is taken 
to be of the form Q oc R~ q . They consider only axisymmet- 
ric perturbations, and this enables them to carry out a local 
stability analysis and to obtain a local dispersion relation 
relating wave frequency u) and wavevector k = (fcji, 0, k z ), 
in cylindrical polar coordinates (R,(j>,z). Their analysis dif- 
fers from previous work in that they take account of both 
compressibility and geometrical terms. They focus on the 
case k z 3> kit in order to study how the MRI modes are 
affected by the strength of the azimuthal field. 

In such a configuration, one might expect to encounter 
two different types of instability: (i) the MRI, which uses the 
magnetic field (here the poloidal component, which is why 
we require B z 7^ and k z 3> kit) to tap the energy in the 
shear, and (ii) buoyancy-driven instability, which makes use 
of the radial structure of the equilibrium configuration to 
tap the (radial) effective gravity. Pessah & Psaltis find three 
regions of instability. Regions I and II (see their Figures 3 
and 5) seem to correspond to the two types of instability 
mentioned above. 

Type I instability occurs at small k z and cuts off at a 
finite value of k z — ksu- This corresponds roughly to modes 
being unstable only when the time for a vertical magnetic 
wave, with wavespeed vaz = Bz/^/Airp, to cross one (verti- 
cal) wavelength (A ~ l/k z ) is longer than the orbital time 
1/fi. This is the standard criterion for instability to MRI, 
and corresponds physically to the situation in which a ver- 
tical field line is unable to straighten itself fast enough to 
counter rotational effects (centrifugal force). Thus instabil- 
ity requires 



k z VAz 5- fi- 



(6) 



Pessah & Psaltis also find that the MRI modes are sta- 
bilized if the azimuthal field is strong enough. They find 
that modes of Type I are unstable only for low values of 
B$ < Bcrit such that 



Here, as before, VA<t> ~ B^/^/inp is the azimuthal Alfven 
speed, vk is the Keplerian azimuthal velocity and c s is 
the sound speed. The regime we are interested in corre- 
sponds to a suprathermal, strongly azimuthal field, so that 
c s <C VA4, <C vk and B z <C B$. This regime differs from 
the situation in which the predominant field is B z in two 
important respects. First, the azimuthal field exerts an ad- 
ditional restoring force for a k z perturbation, which, for a 
given k z , dominates at large B$. Second, the length of a field 
line corresponding to one vertical wavelength (A ~ l/k z ) is 
now 



d ~ A 



'~B~z 



(8) 



The analysis of Pessah & Psaltis appears to suggest that 
the MRI is no longer able to operate when the azimuthal 
field is increased to the extent that the timescale on which 
the restoring force due to the field operates (~ R/va<p) is 
shorter than the timescale on which pressure equilibrium 
can be established along the field line (~ d/c s ). In this case, 
instability requires 



k VA > VM 
KzVAz -1, ~- 

C S R 



(9) 



2 <r 

V A<t> Z 



VrC s 



(7) 



Putting conditions © and ((9J together, we obtain the con- 
dition for instability given by inequality (0 above. 

Pessah & Psaltis suggest that Type II instability corre- 
sponds to buoyancy modes, in line with the results of Kim 
& Ostriker (2000) for the limit c s -> 0. Modes of Type I and 
Type II are clearly linked in some way, as they undergo an 
exchange of stability at the critical azimuthal field strength. 
Type II modes are unstable only when both inequalities (JS| 
and @ are violated, and this only occurs for large values of 
the azimuthal field strength B^, > _B cr it. Thus, for a given 
vertical field strength B z , MRI modes — which tap the shear 
energy — are able to operate for small k z and B^ < Bcrit, 
whereas for B<f, > B C rit, radial buoyancy modes, with large 
k z , take over. 

The derivation of the stability criterion by Pessah & 
Psaltis requires the presence of a vertical field {B z 7^ 0), 
however small. But analogous behavior is found for the 
case in which the field is purely azimuthal (Terquem & 
Papaloizou 1996). Since MRI draws its energy from the 
background shear, it is evident that instability can arise in 
this case only for nonaxisymmetric modes. Terquem & Pa- 
paloizou use trial local displacements to show that the linear 
operator describing the evolution of the linearized equations 
has a dense or continuous spectrum of oscillation frequencies 
u> that satisfy a local dispersion relation. Their analysis is 
restricted to modes that are essentially incompressible in the 
(R, z)-plane[3 This eliminates the fast MHD mode, which in 
any case is not involved in MRI, and means that the disper- 
sion relation is only a quartic in u). 

In the limit we are considering, with k z kn and c s <C 
va <C vk, the Terquem & Papaloizou dispersion relation 
simplifies considerably. For axisymmetric [m = 0) distur- 
bances, u) is real, implying stability. For non-axisymmetric 
disturbances, we consider first modes with < m <vk /va- 
For these modes the timescale for an azimuthal magnetic 



1 Note that this restriction still permits the effects of compress- 
ibility to operate in the azimuthal direction. 
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wave, with wavespeed ua, to cross one (azimuthal) wave- 
length (~ R/m), is longer than the orbital time 1/S1 (cf. 
condition [6])- For this range of m there are two low fre- 
quency roots which can be unstable (in addition to two sta- 
ble roots with ui + mQ ~ ±fl). What is interesting for our 
current discussion is that the criterion for instability corre- 
sponds approximately to condition (0 derived by Pessah & 
Psaltis. The growth rates are of order ~ mc„/R, and the 
(real parts of the) frequencies in the rotating frame are such 
that W + mfl ~ mv\/ R 2 Q, < va/R- 

For higher values of m, such that m^VK /va, one finds 
instability provided that m<iu/c s . Putting these two con- 
ditions together thus implies that for instability we require 
v\ > c s vk, which is the reverse of condition 0. This may 
correspond to the exchange of stabilities found by Pessah 
& Psaltis between their modes of Type I and Type II. The 
growth rates for these high m modes are ~ va/R, and the 
real parts of the frequencies in the rotating frame are small, 
uj + rafi ~ Q/m. The fact that these modes are almost sta- 
tionary in the corotating frame and that the growth rates 
depend on the field strength suggests, as for the Pessah & 
Psaltis Type II modes, that they are driven by buoyancy. 

The analogy between the findings of Pessah & Psaltis 
for axisymmetric modes in the limit of < B z <C B^,, and 
the analysis of Terquem & Papaloizou for nonaxisymmetric 
modes when B z — 0, seems to imply that the same physical 
processes are at work in both cases. 

The basic conjecture of this paper is that amplifica- 
tion of the toroidal magnetic field occurs only while MRI 
is available to drive a dynamo. Since MRI dies out as B C rit 
is approached from below, we assume that B^ saturates at 
roughly this level. Indeed, the growth rates of both Type I 
and Type II modes are small in the vicinity of B C rit, where 
the exchange of stability between MRI and buoyancy modes 
occurs. However, even if B^ managed to cross this bottle- 
neck and buoyancy modes were excited, this does not mean 
that B^, will continue to grow. Since buoyancy modes do not 
tap the shear, they are unlikely to drive a dynamo capable 
of amplifying the toroidal magnetic field. 



4 CONSEQUENCES FOR DISC STRUCTURE 

In this section we develop a one-zone model for vertical 
disc structure, under the assumption that the toroidal field 
reaches its limiting strength according to the Pessah & 
Psaltis (2005) analysis. The limiting field strength is related 
to the ability of sound waves to propagate along toroidal 
field lines. Because of the large radiative diffusivity of ac- 
cretion discs, such waves are unlikely to be mediated by ra- 
diation pressure, even where radiation dominates the total 
pressure. Therefore we assume that the appropriate sound 
speed to use in estimating ha is that due to the gas pres- 
sure, c g ~ (Pg/p) 1 ^ 2 , and we take va ~ (cqVk) 1 ^ 2 - For a 
thin disc we have the ordering c 9 < ua « hk- When va 
exceeds the sound speed associated with the radiation pres- 
sure, c r ~ (p r /p) 1,/2 , the scale height is given by 



H 
R 



VA 
VK 



(-) 



1/2 



(10) 



instead of the Shakura-Sunyaev result, Hss/R ~ c 3 /vk = 
(c 2 , + c 2 ) 1 ^ 2 /v-k (eq. [3j . If we assume an a— model viscosity 



with kinematic viscosity given by v — aHvA (rather than 
the usual v = aHc s , but preserving the assumption that the 
stress is ~ aps with ps being essentially the total pressure), 
then the inflow speed (for a Keplerian disc) is given by 



3 (H\ 2 
Vi »~2 a \R) V 
and the column density is 



3 

2 a ° 9 



M 



(11) 



(12) 



where M is the local mass accretion rate through the disc. 

It is convenient to normalize radii to the gravitational 
radius, R g = GM/c 2 , and the accretion rate to the Edding- 
ton accretion rate, which we define as MEdd = L-^dd/c 2 = 
\-kGM I kc, where k, is the opacity and M is the mass of 
the central object. Defining m = M/M^dd and x = R/R g , 
we obtain an expression for the transverse optical depth 
through the disc, 

4 m c _! 

T = IjK ~ X 



3ac, 



(13) 



If the disc is radiative, the flux from each side (at x 3> x^x, 
where Xi n = 6 is the radius of the innermost stable orbit for 
a Schwarzschild black hole) is given by 



F 



3 GMM 



2p r c 



(14) 



8tt R3 2GMkx 3 r ' 
where p r is the radiation pressure near the midplane and the 
factor 2 in the last relation arises because the flux escaping 
from each side of the disc traverses half the optical depth 
(which is assumed to be > 1). We now solve for the radiation 
pressure, 



(15) 



The density inside the disc is given by 

3/2 



P = 



2H 3 GMk 



—x~ 9/4 . 
a 



(16) 



The results presented so far are fully general. To obtain 
numerical results from them, however, one must determine 
the value of c g and check whether va > c r . This depends on 
the equation of state, the opacity, and the relative impor- 
tance of radiation and gas pressure, which we now address. 

4.1 The inner regions of black hole accretion discs 

4-1.1 The disc structure 

We consider here the inner regions of a standard Shakura- 
Sunyaev disc around a black hole where radiation pressure 
dominates over gas pressure at small radii for values of m 
that are not too small. The opacity is dominated by electron 
scattering and we take k = K es — 0.4 cm 2 g _1 . We begin by 
calculating the magnetically dominated disc structure under 
the assumption pb > p r - We then check the self-consistency 
of this assumption a posteriori, to find the radius within 
which p r > pb and magnetic support can be neglected. 

To determine both c 9 and p r , we assume that the radia- 
tion inside the disc is in local thermal equilibrium (LTE) and 
set T — Tlte = (3p r /a) 1,/4 . We then use eq. (|15fl to solve for 
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T with c g — (fcT//i) 1//2 , where p « 0.6m p is the mean mass 
per particle. We verify the LTE assumption a posteriori. We 
find that i>a > c r for x > 20(am) ' rh ' , implying that 
magnetic pressure dominates over other forms of pressure at 
all radii, except possibly for the region immediately outside 
a black hole or neutron star. Radiation pressure dominates 
over gas pressure for x < 2.0 x 10 3 (am) 2 / 13 m 8 ' /13 , but since 
magnetic pressure determines the vertical structure through 
most of this region, there is no difference in properties be- 
tween the radiation and gas pressure-dominated zones. 

The electron scattering optical depth to the disc mid- 
plane is 

r cs ~ 2.5 x 10 2 cr 8/9 m 1/9 m 7/9 :c- 5/9 . (17) 

To check the assumption of LTE, we also need the absorption 
opacity, which we take to be the Kramers bound-free opacity 
for solar abundances, Kbt = 1-6 x 10 24 pT~ 7 ^ 2 cm 2 g _1 . The 
effective optical depth for thermalization is then 

Tsfc = Tcs fOOL) 1/2 „ 0.29a- n/12 m l/12 m 1/3 x 5 ^, (18) 

indicating for typical values of a that in general, LTE is a 
reasonable approximation. 

Finally, we estimate the disc thickness. Substituting 
into eq. (|10|l . we obtain 

»,„ 1n -2/ \-l/18 ■ 1/9 1/36 2 -1 

— ~ 7.3 x 10 (am) ' m ' x ' cm g . (19) 
R 

The disc aspect ratio, H/R (the "opening angle"), is nearly 
independent of radius and is extremely insensitive to param- 
eters, with a magnitude several times larger than a compa- 
rable Shakura-Sunyaev disc without magnetic support. 

4-1.2 Thermalization and Colour Temperature 

Photospheric colour temperatures typically exceed effective 
temperatures in the inner regions of standard accretion discs 
because the radiation is thermalized at a significant scatter- 
ing optical depth. Colour corrections of luminous accretion 
discs, /coi = T co i/r c ff, usually lie in the range ~ 1.5 — 2 
(Shimura & Takahara 1995; Davis et al. 2005), correspond- 
ing to a scattering optical depth ~ 5 — 20 at the thermal- 
ization layer. When r» > 1, the colour correction is given 
roughly by / co i ~ (k cs /k») 1 '' 8 , where «» is the absorption 
opacity evaluated at the local disc colour temperature and 
at the density of the thermalization layer, p*. 

An accurate calculation of the colour correction requires 
a detailed model for the structure of the disc photosphere, 
which is beyond the scope of our one-zone model for the 
vertical structure. However, it is clear that a magnetically 
supported disc should have a higher colour temperature than 
the equivalent Shakura-Sunyaev model (Blaes et al. 2006). 
Here we make a crude estimate of this effect. 

The absorption opacity at the thermalization layer is 
given by 

p* ( T \ 7/2 

ft* = K b f(p*,Tcol) = Kbf — ^ — , (20) 

P VJcol/ 

where the unsubscripted density and temperature repre- 
sent values at the midplane. We also have T/T co i — 
(T cS /T coloul )(T/T cB ) = Te/Vcoi- Substituting for the tem- 
perature ratio in eq. (|20|) and using eq. lfl8)l. we can solve 
for /coi to obtain 



1/4 / \ -2/9 -25/1 

since p*/p < 1. In the inner portions of the magnetically 
dominated zone, eq. (|21[l predicts a colour correction sub- 
stantially larger than that of a standard disc. 

4-1.3 Thermal and Viscous Stability 

A major problem with standard a-disc models is that they 
are viscously and thermally unstable in their inner regions 
where electron scattering dominates the opacity (Lightman 
& Eardley 1974; Shakura & Sunyaev 1976; Pringle 1976). In 
contrast, the magnetically supported disc models we propose 
here are not subject to these thermal and viscous instabil- 
ities. In all radiative a— disc models in which the opacity 
is dominated by electron scattering, the dissipative heating 
rate scales as Q + oc H 2 Y, and the radiative loss rate scales as 
Q- oc p r /E oc T 4 /E, where we assume LTE. From eq. pop , 
we find T oc H 4 at fixed R, implying that Q~ oc H le /E. At 
constant E, losses increase with H much more rapidly than 
heating, implying that the magnetically supported discs are 
thermally stable. This result can be contrasted to a radi- 
ation pressure-dominated a-disc, in which T oc H 1 ^, im- 
plying instability. (Gas pressure-supported a-discs, on the 
other hand, are thermally stable.) 

Now consider viscous instability. On viscous timescales, 
thermal balance is maintained, implying Q + — Q~ . There- 
fore, E oc H 7 for the magnetically supported disc. The vis- 
cous couple satisfies G oc _ff 2 E oc E 9//7 . Since the viscous 
couple is an increasing function of surface density, the disc 
is stable (Lightman & Eardley 1974; Pringle 1981). 

4.2 The thickness of CV and LMXB discs 

By using the same analysis as above, but now using 
Kramers, rather than electron scattering, opacity, we can 
calculate the thickness of the outer regions of discs in bi- 
nary stars. As a specific example, we consider a disc around 
a 1 Mq white dwarf, at a radius of R — 10 10 7?io cm, accret- 
ing at a rate of M = 10 18 Mi 8 g s _1 w 10" 8 M Q yr _1 . This 
corresponds to the outer regions of the disc of a dwarf nova 
in outburst or of a nova-like variable. Under the standard as- 
sumptions the thickness of such a disc would be H/R ~ 0.05 
(see, for example, Frank, King & Raine 2002). However, if 
we use the assumption of strong toroidal fields put forward 
in this paper we find that the disc thickness is given by 

§ = QMa-^MT ) _15/68 RT- (22) 

From this we see that not only are the strongly mag- 
netized discs proposed here thicker than standard ones by 
about an order of magnitude, but they are also slightly more 
flared, with H/R oc i? 0132 rather than the usual oc R 0125 . 
Thus such discs are more subject to irradiation from a cen- 
tral luminosity source. 

4.3 Reduction of Disc Self-Gravity 

Magnetic pressure support can reduce the effects of self- 
gravity in the outer parts of an accretion disc by thickening 
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the disc and reducing its density (Pariev et al. 2003). For a 
disc supported by gas pressure, local gravitational instability 
and fragmentation are expected to limit the accretion rate 
to 

M < M m ~ = 5 x 10- 4 a2 , 1 3 / 2 M Q yr" 1 , (23) 

where T 100 = T/100 K (Shlosman & Begelman 1987). In a 
galactic nucleus, far from the black hole, the disc tempera- 
ture is set by environmental influences (external irradiation, 
cosmic rays, etc.), rather than the internal dissipation in the 
disc (Shlosman & Begelman 1989). 

The maximum accretion rate given by eq. (|23p is too 
small to power luminous active galactic nuclei (AGN). 
Magnetic pressure support, according to our prescription, 
would increase this upper limit by a factor (vk/c 9 ) 3 ^ 2 ~ 
10 5 (M 9 /Tioo.Rpc) 3/4 , where black hole mass and radius are 
in units of 1O 9 M0 and pc, respectively. In Eddington units, 
the maximum accretion rate that could pass through a ra- 
dius R would be increased to 

m max ~ 25aM 9 - 1/4 T 1 3 / 4 i?- 3/4 . (24) 

The decline of rh ma x with R continues out to a radius Rbh, 
where the gravitational potential of the galaxy is comparable 
to that of the black hole. If the galactic nucleus is represented 
by a singular isothermal sphere with velocity dispersion a, 
the "bottleneck" occurs at Rbh ~ GM/2a 2 — beyond Rbh, 
the carrying capacity of the disc, at fixed T, is independent 
of radius. If we use the M — a relation (Ferrarese & Merritt 
2000; Gebhardt at al. 2000; Tremaine et al. 2002), M 9 ~ 
0.13(<r/200 km s" 1 ) 4 , to eliminate a in favor of M, we obtain 

rh max ~ 3aM~ 5/8 T^o ■ (25) 

Thus, magnetic support may permit accretion at close to the 
Eddington limit, provided that a and Tioo are not too small. 
(Note that we defined MEdd without an efficiency factor, so 
rh ~ 10 is required to produced an Eddington luminosity 
with an efficiency of 0.1.) 



5 DISCUSSION AND CONCLUSIONS 

We are well-aware of the speculative nature of some of the 
assumptions that have gone into our proposed disc model. 
Our proposal that da saturates at ~ (csWk) 1 ^ 2 is an ed- 
ucated guess based on a very simple interpretation of the 
likely conditions under which MRI is able to drive strong 
turbulence. We implicitly assume that MRI operates at all 
heights above the disc midplane, amplifying B^ in situ; how- 
ever, may also be advected to high latitudes by buoyancy. 
If buoyant transport from below dominates, then va could 
be even larger than we suggest. 

The surviving MRI modes when va 2> c s have large 
k z , implying that the resulting turbulence would be driven 
on small scales. This is in contrast to the more common 
assumption in MRI calculations, that va <C c s , in which 
case turbulence is driven on scales up to the disc thickness. 
In order to build up the strong, large-scale toroidal field in 
our picture, the dynamo might have to involve an inverse 
cascade, the possibility of which is by no means certain. If 
the production of a large scale toroidal field is not efficient 



in this limit, then our proposed field strength could be a 
considerable overestimate. 

Alternatively, we have suggested the possibility that ac- 
cretion discs with strong toroidal fields of the kind we discuss 
here may form a second stable branch of accretion disc con- 
figurations, perhaps only accessible through some sudden 
change in properties — for example, a thermal instability 
(cf. Section 3). This kind of bimodality of thin discs could 
help to explain the spectral differences of otherwise similar 
accreting white dwarf systems in CVs (cf. Section 2.1). 

Whether or not our quantitative estimate of the satu- 
rated field strength is accurate and physically realized, we 
wish to emphasize the role that a magnetic "equation of 
state" — a relation between a suprathermal B<f, and other 
fluid variables — could play in altering the structure and 
stability properties of a— discs. 
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